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Abstract 

The peculiarities of thermal expansion of bcc and fee zirconium films with (100) and (110) crystallographic orienta¬ 
tions are studied at a constant zero pressure by the molecular dynamics (MD) method with a many-body interatomic 
interaction potential obtained in the embedded atom model. It is shown that after relaxation the cubic lattices become 
tetragonal {bet and fet), and for the metastable fet films the linear coefficients of thermal expansion in the film plane 
have a negative value in a wide temperatuThe peculiarities of thermal expansion of bcc and fee zirconium films with 
(100) and (110) crystallographic orientations are studied at a constant zero pressure by the molecular dynamics (MD) 
method with a many-body interatomic interaction potential obtained in the embedded atom model. It is shown that 
after relaxation the cubic lattices become tetragonal {bet and fet), and for the metastable fet films the linear coefficients 
of thermal expansion in the film plane have a negative value in a wide temperature range. The total and local vibrational 
density of states (VDOS) polarized along the x, y, z axes is calculated for the surface and interior layers of bet and fet 
Zr films as a function of the temperature. It is shown that the peculiarities of the behavior of the vibrational density of 
states of surface atomic layers manifest themselves in the anisotropy of the changes of the film lattice parameters with 
temperature variation. A decrease in the lattice parameters with increasing temperature is observed in the directions 
where there occurs a softening of the local vibrational modes, re range. The total and local vibrational density of states 
(VDOS) polarized along the x, y, z axes is calculated for the surface and interior layers of bet and fet Zr films as a 
function of the temperature. It is shown that the peculiarities of the behavior of the vibrational density of states of sur¬ 
face atomic layers manifest themselves in the anisotropy of the changes of the film lattice parameters with temperature 
variation. A decrease in the lattice parameters with increasing temperature is observed in the directions where there 
occurs a softening of the local vibrational modes. 
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1. Introduction 

In recent years, an ever increasing number of papers 
have been devoted to the study of so-called “negative” 
materials having negative coefficients of thermal expan¬ 
sion [1, 2, 3, 4]. In 1996 the research team led by John 
Evans published the paper [I] in which a negative co¬ 
efficient of thermal expansion (CTE) was found for the 
ZrW 20 s ceramics in a wide range of temperatures from 
0.3A to 1050A. Later many other compounds with simi¬ 
lar properties, named Negative Thermal Expansion (NTE) 
materials, have been discovered [5]. 

These materials with unusual properties are interesting 
not only theoretically, but also from an industrial stand¬ 
point. Materials with controlled thermal expansion can 
find several practical applications, such as switches and 
sensors, while materials with zero thermal expansion in a 
wide temperature range will be used in high-precision en¬ 
gineering. In instrument-making industry materials with 


Email address: elenaaftiudm.ru (E.B. Dolgusheva and V.Yu. 
Trubitsin) 


strictly regulated values of linear CTE in specified tem¬ 
perature ranges of operation are in great demand. Re¬ 
cently the Evans team has created a material whose prop¬ 
erties can be controlled throughout the thermal expan¬ 
sion range: from positive through zero to negative CTE 
[6]. The composite oxide ZrMo20s is a well studied com¬ 
pound distinguished by a negative CTE, whereas the tin- 
containing analog of this cubic material, SnMo 20 s, ex¬ 
pands when heated. The authors of that paper state also 
to have demonstrated the possibility of introducing var¬ 
ious amounts of zirconium into the lattice so that the 
obtained ceramic material can have a positive, negative 
or zero CTE. In essence the “positive” nature of the tin- 
containing material is balanced by “negative” zirconium, 
both metals being key elements of the same crystalline 
structure. Additionally, in Refs. [7, 8] it was shown by in¬ 
elastic neutron scattering that in NTE structures a soften¬ 
ing is observed in the VDOS curve. Thus a negative ther¬ 
mal expansion coefficient is observed in materials which 
are in a metastable state (for example, near a structural 
transition), have a strongly anharmonic character of in- 
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teratomic interaction, and/or a softening of the phonon 
modes is evidenced in the VDOS curve, which results in a 
negative value of the Griineisen parameter. 

Zirconium is known to be a strongly anharmonic metal. 
As we see, it plays a decisive role in the formation of 
negative GTE in the NTE structures. We believe, then, 
that the study of the thermal expansion of Zr nanofilms 
is an important task, because elucidation of the depen¬ 
dence of linear GTE on the temperature and film thick¬ 
ness will make it possible to reveal the peculiarities of the 
formation of thermodynamical properties of nanofilms and 
draw some conclusions about the influence of the surface 
structure on their properties. Clarification of the nature 
of special physical features of nanosystems is of interest 
for both fundamental problems of the condensed matter 
physics and possible practical applications. 

Earlier [9], we constructed a phase diagram 
Temperature-Thickness for zirconium films. This di¬ 
agram presents the regions of final structures stable for 
no less than 1 ns, one time step being equal to 1 fs. Here 
we should mention some important points. Firstly, the 
size effect was detected. Two characteristic sizes may be 
noted. The first one is of 5 unit cells 2nm): up to this 
size the film just falls apart, and the second critical size is 
the film thickness of 17 unit cells (~ 6.1nm) at which the 
mechanism and sequence of structural transformations 
change: from this size on the bcc (001) film remains 
stable in a sertain temperature range. Secondly, what 
is especially important, structures not observed in pure 
Zr in the usual bulk state were found. These are fee 
and orthorhombic structures. Recall that at atmospheric 
pressure only two phases are observed in Zr bulk samples: 
a high-temperature bcc phase, and an hep phase found 
below 1136Ar. Thus, zirconium provides a possibility 
of studying the properties of a strongly anharmonic 
material existing in both stable and metastable states, in 
conditions close to structural transformations. 

The molecular dynamics method makes it possible to de¬ 
termine the linear GTE in different crystallographic direc¬ 
tions, namely, in the direction normal to the film surface, 
and in the plane of the film atomic layers, which allows 
the study of the him thermal expansion anisotropy. Fur¬ 
thermore, this method enables the investigation of physical 
characteristics, such as the lattice dynamics, separately for 
the surface and interior atomic layers of the him. In this 
work the MD method is used to study the peculiarities of 
thermal expansion of the zirconium thin hlms with cubic 
{bcc, fee) crystal structures for different crystallographic 
orientations of the surface at a constant zero pressure. We 
also performed a comparison of the temperature behavior 
of the lattice parameters and vibrational density of states 
of surface and interior atomic layers, polarized both in the 
him plane and in the direction normal to the him surface. 
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Figure 1: Vibrational density of states of bcc Zr. (a) calculated 
for a crystallite with cyclic boundary conditions at T = 900 K; (b) 
experimental values [14]. 


2. Calculation method 

The physical characteristics of zirconium hlms were 
studied by the MD method using the standard XMD pack¬ 
age [10]. The many-body potential (#2) from Ref. [11] con¬ 
structed within the embedded atom model [12] was cho¬ 
sen to describe the interatomic interaction in zirconium. 
In Ref. [11] it was shown that this potential allows one to 
obtain, to a high degree of accuracy, the bcc and hep lat¬ 
tice parameters of zirconium, cohesive energy, elastic con¬ 
stants, melting temperature, and other physical character¬ 
istics. The bulk phonon dispersion curves calculated with 
this potential along the symmetrical directions of the Bril- 
louin zone of bcc zirconium at different temperatures were 
obtained in our paper [13]. A comparison with the experi¬ 
mental data shows that the potential chosen allows one to 
reproduce the experimentally observed features of the Zr 
phonon spectrum including the softening of the transverse 
N-phonon with decreasing temperature, and so it may be 
successfully used in calculating vibrational and thermal 
characteristics of zirconium films. 

In this paper we present only the vibrational density of 
states for bulk zirconium (see Fig.l) as compared with the 
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experiment [14]. The total VDOS, g(w), was calculated 
as the Fourier transform of the autocorrelation velocity 
function, the velocity was averaged over all atoms on an 
interval of 10 ps after a relaxation period of 100 ps at a 
given temperature: 


9{w) 


' E^i(vi(0)|vi(0)) 


exp(iwt). 


( 1 ) 


Here Vi(t) and Vi(0) are the velocities of the i-th atom at 
time t and the initial time, respectively. < > denotes the 
averaging over time and the sum is taken over all atoms. 
For the calculation of the local VDOS (for example, the 
surface atoms VDOS) in the sum were considered only 
the velocities of the surface layer atoms. The local VDOS 
then were multiplied by the weight coefficient. To obtain 
the polarized VDOS all the velocities were multiplied by 
the polarization vector p in calculating the autocorrelation 
function: 


N 

i=l 

As seen from Fig.l, the positions of the main peaks are 
in good agreement with the experiment. 

In calculating the lattice parameters, the coordinates of 
all the crystallite atoms were first averaged over a time in¬ 
terval of 50 ps, then the average values were taken over the 
atomic layer, thereafter the distance between the atomic 
layers in the given direction (x, y, z) was found. 

Note that in calculating the interlayer distance along the 
z axis the near-surface layers (5 layers below the surface on 
each side) were excluded from consideration because of the 
fact that after relaxation near the film surface significant 
deviations from the interlayer distance values in the film 
interior are observed. 

The coefficients of thermal expansion were calculated as 
follows: 




( 3 ) 


where /3 is the film GTE; V is the average volume per 
atom only in the film interior at the initial temperature; 
dV is the relative change of this volume on heating the 
crystallite by dT degrees. For the calculation of the linear 
GTE from our simulation we used the following expression: 




1 dLx,Y,z 
Lx,y,z dT 


( 4 ) 


Here is the linear GTE in the direction normal to 
the film plane along the 2 axis; ax,y are the linear GTE 
in the film plane along the x and y axes; L is the average 
distance between the atomic layers at the initial tempera¬ 
ture along the direction chosen; dL is the relative change 
of this distance on heating the crystallite by dT degrees. 


3. Lattice Parameters and Vibrational Density of 
States 

Note at once that the lattice cubic symmetry in bcc and 
fee films changes after relaxation to tetragonal. As a re¬ 
sult, the bee films turn into bet, and fee into fet ones. 
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Figure 2: Lattice parameters a and c for the films after relaxation: 
(a) for fee films: (b) for bee films. The experimental values for fee 
Zr film were obtained: in [15, 16] (exp-1); in [17] (exp-2), and in [18] 
for a bulk bee Zr sample. 


Figure 2 presents the lattice parameters a and c aver¬ 
aged over 50 ps for films of different thickness indicated 
on the abscissa in unit cells (u.c.) of the corresponding 
structures. For fet films the calculation was performed 
at r = SOOAT (Fig.2(a)), and for bet films at T = OOOAT 
(Fig.2(b)). The statistical error was 0.001 A at SOOAT and 
0.002 A at 900iF. The lattice parameters for fet films with 
(001) and (110) crystallographic orientations are rather 
close in magnitude, while for bet (001) and (110) films 
they differ significantly. The bet (001) film is stabilized 
only from a thickness of 17 unit cells {N^ = 17m. c.). In 
the same figure are shown the experimental lattice param¬ 
eters obtained in Refs. [15, 16] (4.6lA indicated by black 
square), and Ref.[17] (4.63A indicated by black point) for 
fee films, and in Ref. [18] (3.609A) for a bulk bcc zirco¬ 
nium structure (fig.2(b)). The tetragonal distortion pa¬ 
rameter, c/a, decreases with increasing film thickness. For 
fet films it varies from 1.042 to 1.011. The tetragonality 
of fet zirconium films was also observed experimentally 
[16] with a tetragonal distortion parameter c/a equal to 
1.0220 ±0.022. 

The thickness has a determining effect on various prop¬ 
erties of the films, which is due to their structural peculiar¬ 
ities and the arising mechanical stresses. We calculated the 
averaged distances between atomic layers in films of differ¬ 
ent thickness along the direction normal to the film plane, 
i.e. from one free surface to the other. The size of the ba¬ 
sic crystallite along the x and y axes in all cases was equal 
to 24 u.c. The bet films were considered at T = 900iF, and 
the fet ones at T = 300AT. 
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As shown by the calculations, significant variations in 
the average distance between atomic layers are observed 
in the interface regions of all films after relaxation. These 
areas with strong distortions were eliminated from further 
calculations of the lattice parameters and, accordingly, of 
the coefficients of thermal expansion. In all the cases con¬ 
sidered, five interlayer intervals were assigned to the in¬ 
terface region. In calculating the lattice parameters un¬ 
der changes of temperature, only the values obtained for 
the film interiors were taken into consideration. In this 
way the lattice parameters of bet and fet films with crys¬ 
tallographic orientations (001) and (IIO) were calculated 
at different temperatures. At the same time, we calcu¬ 
lated and analyzed the distribution of the total and local 
vibrational densities of states and their temperature de¬ 
pendences in order to determine the extent to which the 
structural stability of bet and fet Zr films is affected by 
the interior and surface vibrations. To this end the sur¬ 
face atomic layers were separated out and all calculations 
were performed separately for surface and interior atomic 
layers. 
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Figure 3: Vibrational density of states: for bet (001) and (110) 
films at T = 1200K (top panels); for fet (001) and (110) films 
at T = 300iC (bottom panels). All the films are 20 u.c. thick. 
The contribution of vibrational states from the surface atomic lay¬ 
ers (multiplied by 2) is indicated by the red line, from the interior 
atomic layers by the green line; the black line shows the vibrational 
states of all the atoms of the film. 


The top panels of Fig.3 present the vibrational density 
of states for bet films with (001) and (110) surface orien¬ 
tation, and the bottom panels for fet films with the same 
surface orientation. As can be seen, for the bet (001) film 
all surface vibrations marked by the red line are in the 
low-frequency region, which is indicative of the dynamic 
instability of the surface, while for the film with (110) sur¬ 
face these vibrations are uniformly blurred throughout the 
frequency range pointing to the dynamic stability of the 
surface. In both cases VDOS for the whole film are prac¬ 
tically identical. 

Returning to the phase diagram Temperature-Thickness 


[9], we can assume that the critical size of 17 u.c (6.1nm) is 
due to the fact that with an increase of the film thickness 
the proportion of surface atoms decreases and so does the 
contribution from the surface vibration modes. As a result 
the film with (001) surface becomes energetically more fa¬ 
vorable, whereas calculations of the bee lattice energy at 
T — OAT showed the (110) film to be energetically favorable 
for all the film thicknesses (see Fig.6 in Ref.[9]). 

In the /ct(OOl) film surface vibrations give a small con¬ 
tribution to the high-frequency region too, while in the 
/ct(llO) film all surface vibration modes are concentrated 
in the low-energy region. Our calculations on heating the 
films show the /ct(OOl) film to be stable over a wider 
temperature range: it undergoes structural transformation 
only at lOSOAT, whereas the /ct(llO) film even at 900Ar. 
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Figure 4: Local VDOS for the surface layers of the /ct(OOl) Zr film 
polarized along the x, 2 : axes at different temperatures. 


Figure 4 shows the results of calculation of the local 
polarized VDOS for surface atoms of a /ct(OOl) film of 
thickness Nz = 20u.e. as a function of the temperature 
indicated on the right. Here one can see noticeable shifts 
of the low-frequency peaks from the energy value of 7meV 
to 2.5meV for local VDOS along both x and y axes. At 
the same time the contribution of high-frequency (about 
ISmeV) vibrations polarized along these axes substan¬ 
tially decreases. The position of the peak corresponding to 
the surface atom displacements along the z axis on heat¬ 
ing first shifts slightly to the region of higher energies, and 
starting from the temperature T = 700K this peak shifts 
in the opposite direction. This behavior, in our opinion, 
is due to changes in interlayer distances in the interface 
region. At low temperatures oscillations of the interlayer 
distances are observed, whereas with increasing tempera¬ 
ture these oscillations disappear, and all distances between 
the atomic layers in the interface region (5 atomic layers 
from the surface) increase. Just such a reconstruction oc¬ 
curs around 700K. Before the structural transformation 
(at T = lOOOAT) all the atomic layers in the near-surface 
area are spaced apart at larger intervals than the layers of 
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the film interior. 
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Figure 5: Local VDOS for the surface layers of the /ct(llO) Zr film 
polarized along the x,y, z axes at different temperatures. 

The results of the local polarized VDOS calculation for 
surface atoms of a fct (110) film are presented in Fig.5: 
here the shift of the low-frequency peaks is clearly observed 
only in one direction along the x axis. Note that different 
nature of the atomic layer displacement in the fct (110) 
and (001) films leads to the formation of different struc¬ 
tures. The fct (110) film transforms into an hep (0001) 
at 900K, and the fct (001) film into a bcc (110) one at 
a temperature of 1050iF. Note also that on heating the 
oscillations of interlayer distances persist in the fct (110) 
film up to the temperature of structural transition. 

We also calculated the temperature dependences of the 
lattice parameters ax,ay, and c of fct films 20 u.c. thick 
(see Fig.6). Red color shows the changes in the /ct(OOl) 
lattice parameters. One can see that the parameter c 
increases with temperature, while ax and ay decrease 
equally in the whole temperature range considered. In the 
/cf(llO) film (blue color) only the parameter Ux decreases, 
whereas ay and c grow. 

Also were calculated the temperature dependences of 
the local VDOS polarized along the x,y, and z axes for 
both all interior atomic layers of fct films and only two 
interior layers. However, no shift of vibrational modes to 
the low-energy region was detected in either case. Thus, 
for metastable films with fct structure there is a clear cor¬ 
relation between the temperature behavior of the surface 
atom vibrations and the lattice parameters. 

Analogous VDOS calculations were performed for bet 
zirconium films of thickness = 20u.c. with (001) and 
(110) surface orientation. In the bct{110) film, on heat¬ 
ing, all the lattice parameters ax,ay,c increase, and no 
softening of vibrational modes is observed in the tempera¬ 
ture behavior of the local phonon densities polarized along 
the X, y, z directions for both surface and interior atomic 
layers. In the &cf(001) film, when heated from 600K to 
1200Ar, all the lattice parameters also increase, but on 



Temperature, K 

Figure 6: Temperature dependence of the lattice parameters of fct Zr 
films of thickness 20u.c.: red line for /c£(001); blue line for /c£(110). 


further heating to 1500AT (i.e. as the boundary of the 
structural reorientation transition is approached) a slight 
decrease in the ax and ay parameters is observed, though 
the parameter c continues growing. As seen from the 
phase diagram Temperature-Thickness [9], the bet (001) 
film of thickness = 20u.c. is stable in the temperature 
range from 600K to ISOOAT, and at temperatures above 
ISOOAl it undergoes an orientational transformation into 
bet (110). The decrease in the Ux, ay lattice parameters of 
the 6cf(001) film in the range 1200 — ISOOAT may be asso¬ 
ciated with stresses that arise in these directions resulting 
in a structural transformation. 

4. Coefficients of thermal expansion 

As mentioned above, in the bet (001) film the parame¬ 
ters Qx and ay first increase and then (from 12001^ on) de¬ 
crease. The appearance of such two-dimensional stresses in 
the bet (001) film with an increase in temperature results 
in structural instability. It is obvious that such changes 
in the lattice parameters result in nonlinear variation the 
coefficients of thermal expansion for films, as opposed to 
bulk samples in which these coefficients vary linearly. In 
Fig. 7(a) are presented the CTE for the zirconium films bet 
(001) (red) and bet (110) (green). The figure also shows 
the CTE calculated for a bcc crystallite with cyclic bound¬ 
ary conditions (grey) and the experimental data for a bulk 
Zr sample from Ref. [19] (black). As seen from the fig¬ 
ure, the calculated and experimental CTE values for the 
“bulk” case are in good agreement. 

The linear coefficients of thermal expansion for bet (110) 
and (001) films are presented in Fig.7: along the direc¬ 

tion normal to the film plane (b); Oix,y for the bet (110) 
film (c); the same coefficients for the bet (001) film (d). 

The temperature dependences of the linear CTE were 
also calculated for the fct zirconium films. In Fig.8 the 
linear CTE are shown for the films fct (001) (red) and fct 
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Figure 7: Coefficients of thermal expansion of Zr films as a function 
of the temperature: (a) for bet films (001) and (110) (red and green, 
respectively), for a crystallite with cyclic boundary conditions (grey), 
and the experimental data for the bulk sample[19] (black); (b) linear 
CTE along the 2 : axis for the bet (110)(red) and (110)(green) films. 
Linear CTE along the x and y axes for the films bet (110) (c), and 
bet (001) (d). 


dependences of CTE are presented in Fig.8(d) for the fet 
films (001) (red) and (110) (blue). 

Thus, calculations of the coefficients of thermal expan¬ 
sion and linear CTE for zirconium films from the molecular 
dynamics simulation show that for the stable bet (110) film 
all CTE are positive throughout the temperature interval 
investigated. The linear CTE of the bct{001) film in the 
film plane are alternating-sign, the sign being changed to 
negative just before the structural transformation. 

For metastable fet (001) and (110) films CTE remain 
positive, but their behavior depends substantially on the 
film surface orientation. At high temperature the coeffi¬ 
cient is of essentially nonlinear character for the fet (001) 
film, and also for the fet (110) film the CTE curve exhibits 
a hump. In the whole temperature range the linear CTE 
ax and ay are negative for the fet (001) film, whereas 
for fet (110) only the coefficient ax is negative, and ay is 
greater than zero. 


5. Conclusion 


(110) (blue): along the normal to the film plane (a); 

ax and ay in the film plane(b). 
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Figure 8: Temperature dependence of linear CTE for the films fet 
(001) (red) and fet (110) (blue): along the normal to the film plane 
(a); in the film plane (b). Temperature dependence of the volume 
per atom of fet films (c); temperature dependence of the CTE for 
fet films (d). 


The fet (001) film has negative coefficients of linear ex¬ 
pansion along the axes parallel to the film plane over the 
whole temperature range considered. In the fet (110) film 
the linear CTE is negative only along the x axis, and along 
the y direction it is positive. Note that the fet (110) film 
becomes unstable even at T = OOOAT, there occurs a transi¬ 
tion into an hep structure. The change of volume per atom 
in the film interior is shown in Fig.8(c) as a function of the 
temperature. In both films the volume increases owing to 
a more rapid growth of the parameter c. The temperature 


The calculation of the total and local vibrational densi¬ 
ties of states polarized along the axes x, y, z as a function 
of temperature for the surface and interior layers of bet 
(stable) and fet (metastable) zirconium films has shown 
a correlation between the VDOS of surface atomic layers 
and the anisotropy of the lattice parameters variation with 
temperature. For example, in (100) and (110) fc± Zr films 
the decrease of the lattice parameters with increasing tem¬ 
perature is due to the “softening” of low-frequency peaks 
in the VDOS curves. Furthermore, it is shown that at 
temperatures close to that of the structural transition the 
zirconium films can also exhibit negative linear coefficients 
of thermal expansion. 

It is known that the negative CTE observed in bulk sam¬ 
ples of composite materials (ceramics, alloys) either are 
associated with vibrations of bridging oxygen (Me-O-Me) 
[1], or they arise in highly anisotropic layered structures 
as a result of the membrane effect [20], In the case of zir¬ 
conium thin films neither mechanism of the appearance of 
negative linear CTE is feasible, firstly, because of the lack 
of light atoms necessary for the bridging mechanism, and 
secondly, the existing anisotropy of interaction between 
the atoms in the film plane and in the direction normal 
to the plane does not suffice for the membrane effect to 
arise. Thus, the reasons for the appearance of negative 
coefficients of linear thermal expansion in zirconium films 
require further investigation. 
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